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1.  Introduction 


There  is  substantial  interest  in  block  copolymers  due  to  its  tailorability,  where  morphology  can 
be  easily  manipulated  through  four  factors:  monomer  type,  polymer  architecture,  system 
composition,  and  molecular  size.  Important  quantities,  such  as  mechanical  and  physical 
properties,  are  influenced  by  the  phase  morphology  and  are  essential  when  determining  potential 
applications  (1 ).  Due  to  the  large  parameter  space  associated  with  block  copolymers, 
computational  modeling  can  be  used  to  determine  the  relevant  parameters  associated  with  the 
desired  phase  behavior  in  order  to  speed  up  the  material  design  process. 

One  system  of  interest  is  in  thermoplastic  elastomer  (TPE)  gels  composed  of  triblock  or 
multiblock  copolymers  immersed  in  a  midblock-selective  solvent.  The  copolymer  is  composed 
of  a  rubbery  midblock  and  short  glassy  end  blocks,  where  the  latter  can  cluster  together  to  fonn 
microdomains,  which  act  as  physical  crosslink  sites  situated  in  a  continuous  matrix  composed  of 
the  midblock  (2).  Some  examples  of  TPEs  are  poly[styrene-butadiene-styrene]  (SBS), 
poly[styrene-isoprene-styrene]  (SIS),  poly[styrene-isobutylene-styrene]  (SIBS),  and 
poly[styrene-(ethylene-co-butylene)-styrene]  (SEBS)  (3).  Out  of  these  three,  there  is  particular 
interest  in  SEBS,  which  is  the  hardest  and  most  resilient  against  degradation,  and  is  used  in  a 
variety  of  applications  such  as  biomimetic  gels,  adhesives,  sealants,  coatings,  shoe  soles,  car 
parts,  and  wire  insulation  (3). 

While  computational  modeling  is  ideal  for  investigating  the  block-copolymer  morphology, 
minimal  research  has  been  performed  on  TPE  gels.  Computationally,  this  gel  can  be  modeled  as 
an  ABA  triblock  immersed  in  a  B-selective  solvent.  In  this  framework,  Kim  et  al.  ( 4 )  used  grand 
canonical  Monte  Carlo  to  study  the  shape,  size,  distribution,  and  internal  structures  of  micelles 
formed  at  low  polymer  volume  fractions.  Mechanical  behavior  of  the  micellar  phase  was  also 
studied  by  Sliozberg  et  al.  (5)  using  dissipative  particle  dynamics.  For  a  highly  solvated  system 
(80-90  vol.%)  Shull  and  coworkers  ( 6)  used  the  unit  cell  approximation  in  self-consistent  field 
theory  to  study  the  effect  of  the  triblock  copolymer  midblock  length  and  temperature  on  the 
micellar  size,  geometry,  aggregation  number,  bridging  fraction,  solvent  osmotic  pressure,  and 
form  factor.  Although  insight  into  TPE  gels  has  been  gained  through  these  studies,  they  are 
limited  to  systems  with  high  solvent  concentrations.  For  a  slightly  different  system  of 
polypropylene  oxide-ethylene  oxide-propylene  oxide]  triblock  in  aqueous  solution,  where  the 
solvent  has  a  slight  attraction  towards  the  midblock,  Fraaijie  and  coworkers  (7)  presented  select 
morphologies  for  a  few  polymer-solvent  compositions.  Even  so,  we  are  interested  in  a  more 
comprehensive  study  for  systems  where  the  solvent  has  a  strong  selectivity  towards  the 
midblock. 
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We  have  used  the  dynamic  density  functional  theory  (DDFT)  to  study  the  phase  behavior  of 
SEBS  and  the  midblock  selective  hydrocarbon  oil  tetradecane  (C14)  as  a  function  of  solvent 
volume  fraction  and  the  S-EB  interaction  parameter,  which  is  inversely  proportional  to  the 
temperature.  Additional  calculations  for  the  spherical  morphology  (at  50/50  vol.%  SEBS/C14) 
were  performed  to  detennine  the  effect  of  the  solvent  molecular  weight  on  the  phase 
morphology.  Lastly,  using  a  finite-element  based  mechanical  homogenization  method  for  linear 
elasticity,  we  have  also  calculated  the  Young’s  and  shear  moduli  for  select  morphologies. 


2.  Models  and  Methods 


A  polymer  gel  exhibits  both  thennal  and  environmental  stability,  and  possesses  enhanced 
mechanical  properties  due  to  its  ability  to  form  physical  crosslinks  where  the  triblock 
copolymers  bridges  two  microdomains.  In  this  particular  study,  the  phase  behavior  of  SEBS  in 
C 14  is  studied  using  DDFT,  which  was  developed  by  Fraaije  et  al.  ( 8 ,  9).  The  input  parameters 
required  by  the  method  are  the  number  of  beads  and  the  polymer  interaction  parameters,  which 
can  be  related  to  the  solubility  parameters.  The  SEBS  of  interest  is  manufactured  by  Kraton 
Polymers  under  the  label  Kraton  G1652M,  and  is  composed  of  29.0  to  30.8  wt.%  styrene  with  a 
molecular  weight  of  about  99,000  g/mol.  Although  numerous  methods  are  available  to  calculate 
solubility  parameters,  we  chose  to  use  group  additive  methods  through  the  SYNTHIA  (10) 
module  of  Material  Studios  (MS),  which  uses  the  Fedors  predictions  of  the  solubility  parameters. 
Both  of  the  rubbery  components  and  the  midblock  selective  C14  exhibit  similar  values  of 
17.53  V(J/cm3),  17.1  V(J/cm3),  and  17.53  V(J/cm3)  for  ethylene,  butylene,  and  tetradecane, 
respectively,  which  indicate  that  they  are  miscible  within  each  other  relative  to  the  styrene  block 
(20.1  V(J/cm3)).  Thus,  the  ethylene-co-butylene  midblock  and  the  C14  solvent  can  be  represented 
as  the  same  molecule,  where  the  SEBS  and  C14  blend  are  modeled  as  A3B15A3  +  Bi.  In  this 
discretization,  the  styrene  block  contributes  about  28%  of  the  SEBS  volume  and  assuming  an 
average  density,  molecular  volume,  and  molecular  mass  of  0.92  g/cnr,  48.9  cm  / mol,  and 
42.08  g/mol,  respectively,  for  ethylene-butylene  and  1.05  g/cm  ,  97.0  cm  /mol,  and  104.2  g/mol, 
respectively,  for  styrene,  each  bead  is  calculated  to  occupy  a  volume  of  approximately 
5100  cm  /mol,  where  each  styrene  and  ethylene-butylene  bead  contains  53  and  105  monomers, 
respectively.  The  solvent  is  also  represented  using  the  same  bead  volume,  which  corresponds  to 
23  molecules  of  C14. 

The  MESODYN  module  through  the  commercial  software  MS  (10)  was  used  to  develop  a  phase 
diagram  for  SEBS  +  C14  for  different  volume  fractions  of  the  solvent,  <D>S  and  the  A-B  interaction 
parameters,  Sab  or  Xab-  MESODYN  uses  an  interaction  parameter  Sab  to  characterize  the 
repulsion  between  the  A  and  B  species,  which  is  inputted  into  the  program  as  the  quantity  Sab/v, 
where  v  is  the  mean  field  excluded  volume  parameter.  This  interaction  parameter  is  related  to  the 
Flory  Huggins  parameter  Xab  through  Xab  =  [2sab  -  £aa  -  SBB]/(2vkT)  (8,  11).  In  addition  to  the 


2 


composition  of  the  triblock  and  the  interaction  parameter,  the  degree  of  polymerization,  N,  is 
also  an  important  parameter  since  the  entropy  of  mixing  is  dependent  on  the  size  of  the  polymer. 
To  capture  this  contribution  an  effective  repulsion  parameter,  XabN,  is  defined,  which  is 
calculated  for  our  system  through  XabN  ~  7.82  Sab  /v  mol/kJ,  where  we  consider  a  temperature  of 
T=323K  and  zero  self-interaction  parameters  Saa  =  £bb  =  0. 

The  simulations  were  performed  on  a  periodic  cubic  simulation  box  with  32x32x32  grid  points 
using  the  recommended  MESODYN  parameters  for  the  grid  spacing,  noise,  compressibility 
parameter,  and  Crank-Nicholson  tolerance.  The  simulation  box  was  sheared  normal  to  the  y- 
direction  at  a  rate  of  5.000  e-3,  5.000e-2,  and  1.000e-2,  for  9,000  time  steps  after  5,000  time 
steps  of  equilibrium,  followed  by  12,000  time  steps  of  equilibration,  where  each  time  step  is 
50  ns.  Shearing  was  perfonned  to  speed  up  phase  separation  and  as  a  mechanism  to  escape  from 
metastable  configurations,  where  the  free  energy  of  each  phase  before  (at  5,000  time  steps)  and 
after  (at  26,000  time  steps)  shearing  are  compared  for  the  three  shear  rates.  The  phase  with  the 
lowest  free  energy  was  then  considered  the  equilibrium  morphology.  This  is  similar  to  the 
method  used  by  Zvelindovsky  for  the  spherical  morphology  for  an  AB  polymer  melt  and  ABA 
block  copolymer  and  solvent  mixture  (12).  Although  shear  transforms  the  phase  morphology  into 
the  cylindrical  phase,  the  cylinders  can  break  or  deform  to  different  morphologies  such  as  the 
spherical  or  gyroid  phase  during  the  equilibrium  period  following  the  shearing  process.  This 
reordering  reduces  the  packing  frustration  associated  with  being  constrained  to  a  highly 
energetically  unfavorable  configuration.  Simulations  were  performed  for  <DS  ranging  from  0  to 
0.9  in  increments  of  0.1,  where  Sab  /v  was  varied  from  4  kJ/mol  to  14  kJ/mol  in  increments  in 
0.2  kJ/mol. 


3.  Results  and  Discussion 


The  phase  diagram  produced  using  DDFT  for  SEBS  +  Cm  system  is  shown  in  figure  1,  where  we 
observed  three  classical  microphases:  gyroid,  cylindrical,  and  spherical.  Although  the  phase 
diagram  was  produced  to  describe  our  system  of  interest,  it  can  be  used  as  a  universal  phase 
diagram  for  any  A3B15A3  +  Bi  system.  The  interaction  parameter  can  be  calculated  for  an 
experimental  system  using  the  equation,  Xab  =  Vref{5A  -  5B]2/(RT)  (13,  14),  where  5A  and  5B  is 
the  solubility  parameters  of  species  A  and  B,  respectively,  and  Vref  is  a  reference  molar  volume. 

Morphologies  obtained  using  MESODYN  are  typically  defective  and  metastable,  where  defect- 
free  configurations  are  very  difficult  to  obtain  (9).  Although  shearing  can  be  used  to  escape  from 
metastable  morphologies,  obtaining  ordered  structures  is  still  challenging.  Nevertheless,  the 
morphologies  can  still  be  differentiated  into  three  microstructures:  spheres,  cylinders,  and  gyroid. 
In  our  simulations,  it  is  difficult  to  discern  between  spheres  arranged  in  BCC  symmetry  and  a 
liquid  micellar  phase,  which  are  stable  at  low  and  high  solvent  dilution  (15),  respectively, 
therefore  they  will  not  be  differentiated  in  this  study. 
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Figure  1.  Phase  diagram  produced  using  DDFT  for  an  ABA  triblock  and  B-selective  solvent.  The  lines  are  drawn 
to  guide  the  eyes  and  are  not  actual  phase  boundaries.  G,  C,  S,  and  D  represent  the  gyroid,  cylindrical, 
spherical,  and  disorder  phases,  respectively. 

For  the  neat  ABA  triblock  copolymer  system  (<DS  =  0),  the  phase  sequence  (disorder  —*■  spherical 
—>  cylindrical  — >  gyroid)  is  observed  with  an  increasing  interaction  parameter  or  decreasing 
temperature.  This  is  the  same  sequence  that  was  obtained  using  self-consistent  field  theory  for  a 
symmetric  (A  blocks  equal  in  length)  triblock  by  Matsen  and  Thompson  ( 16)  at  approximately 
the  same  A  volume  fraction.  At  this  volume  fraction,  they  observed  a  XabN  ~13  associated  with 
the  order-disorder  transition  (ODT)  for  a  triblock  copolymer  of  length  2N(16).  We  observe  a 
similar  value,  XabN  ~17,  after  the  doubling  of  the  triblock  copolymer  N  is  accounted  for.  This 
deviation  along  with  shifts  in  the  phase  boundaries  were  also  observed  by  Horvat  et  al.  (if);  they 
were  attributed  to  the  small  number  of  beads  used  to  resolve  the  chain  and  the  nonlocality  of  the 
non-ideal  interactions.  In  addition,  the  DDFT  method  includes  a  random  noise  term  in  the 
equations  associated  with  the  density  evolution,  which  the  self-consistent  field  theory  (SCFT) 
lacks.  Although  the  noise  is  not  capable  of  extending  the  theory  beyond  the  mean  field 
approximation  ( 1 7),  it  may  shift  the  ODT  to  slightly  higher  values  of  XabN. 

Several  experimental  and  theoretical  studies  have  shown  that  the  thermodynamic  behavior  of  a 
symmetric  ABA  triblock  is  very  similar  to  that  of  a  corresponding  AB  diblock  system  obtained 
by  snipping  the  triblock  in  half.  This  is  because  the  entropy  gain  through  this  snipping  is 
relatively  small,  which  allows  the  phase  behavior  for  the  neat  triblock  to  be  governed  by  the 
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same  principles  (16,  18-20).  Due  to  the  asymmetric  of  the  snipped  triblock  copolymer  (A3B7), 
competition  between  the  A  and  B  stretching  energies  stabilize  morphologies  with  curved 
interfaces  and  the  lamellar  phase  is  energetically  unfavorable  as  the  A  volume  fraction  deviates 
from  0.5  {16,  12,  22).  The  sequence  of  the  curved  morphologies  observed  for  our  system  is  due 
to  a  delicate  balance  between  interfacial  tension  and  chain  stretching,  where  the  former  and  latter 
contributions  dominate  at  high  and  low  XabN,  respectively.  At  low  segregation  or  at  high 
temperatures,  the  disorder  phase  is  stable  since  the  repulsion  between  the  A  and  B  segments  is 
weak  and  entropy  dominates.  Upon  cooling  a  disordered  phase,  a  first  order  transition  to  the 
spherical  morphology  occurs.  Due  to  the  asymmetry  exhibited  by  A3B7,  the  system  favors 
morphologies  with  a  higher  interfacial  curvature,  such  as  the  spherical  phase  at  high 
temperatures,  where  the  majority  component  is  allowed  to  relax  at  the  expense  of  the  minority 
block  (23).  This  order-order  transition  was  first  predicted  by  Leibler  {24)  for  an  AB  diblock 
copolymer  using  Landau  theory  and  since  then  has  been  studied  theoretically  and  experimentally 
(25)  for  numerous  block  copolymer  systems.  An  order-order  transition  to  the  cylindrical  phase 
occurs  when  XabN  is  further  increased  and  the  enthalpic  contribution  becomes  more  significant. 
In  this  situation,  the  spheres  deform  and  coalescence  together  to  form  cylinders  arranged  in  a 
hexagonal  lattice  due  to  thennodynamic  instability  of  the  spherical  interface,  where  the  reduction 
in  interfacial  area  lowers  the  enthalpic  contribution  to  the  free  energy  {26,  27).  When  the 
enthalpic  contribution  begins  to  dominate  at  large  XabN,  the  gyroid  phase  becomes  stable  where 
the  geometry  is  capable  of  producing  interfaces  of  uniform  curvature  while  retaining  domains  of 
uniform  thickness  {16). 

Adding  solvent  into  the  system  not  only  stabilizes  the  cylindrical  and/or  spherical  morphology, 
but  also  delays  microphase  separation  where  the  value  of  XabN  associated  with  the  ODT  is 
shifted  to  higher  values.  As  seen  in  figure  1,  the  value  of  XabN  at  which  the  disorder  phase 
transforms  into  the  spherical  morphology  increases  with  increasing  solvent  concentration  and  is 
due  to  the  large  entropy  of  mixing  gained  by  distributing  the  solvent  molecules  randomly  within 
the  system  {28).  For  very  low  volume  fractions  of  solvent,  <1>S  =  0.1,  the  same  phase  sequence 
(disorder  — »•  spherical  — >  cylindrical  — »•  gyroid)  as  seen  in  the  neat  triblock  is  observed,  where 
the  concentration  of  the  solvent  is  not  large  enough  to  destabilize  this  progression.  Further 
increasing  the  solvent  concentration  from  ®s  ~  0.2  to  ®s  ~  0.4,  destabilizes  the  gyroid 
morphology,  while  contracting  and  enlarging  the  window  of  stability  for  the  cylindrical  and 
spherical  phase,  respectively,  with  increasing  solvent  concentration.  For  ®s  ~  0.5  to  ®s  ~  0.8, 
only  the  spherical  morphology  is  observed,  where  microphase  separation  does  not  occur  at 
higher  solvent  volume  fractions.  The  phase  sequence  observed  at  ®s  ~  0.5,  spherical  — >  disorder, 
with  increasing  temperature  was  also  verified  experimentally  by  Sugimoto  et  al.  (29). 

This  phase  behavior  observed  in  figure  1  can  be  explained  using  the  trajectory  approach 
described  by  Hanley  et  al.  (30)  and  Suo  et  al.  (30,  31).  With  increasing  solvent  concentration,  we 
observe  a  phase  transition  (horizontal  trajectory)  from  the  gyroid  — >  cylindrical  — >  spherical  — » 
disorder  at  high  A-B  repulsion  (16),  which  is  the  same  phase  sequence  that  is  observed  for  a  neat 
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triblock  copolymer  when  the  B  volume  fraction  is  increased  from  ~70  vol.%.  This  is  due  to  the 
B-selective  solvent,  which  swells  the  B  domain  and  modifies  the  relative  concentration  of  the  A 
and  B  blocks  such  that  the  effective  B  volume  fraction  increases  compared  to  the  pure  triblock 
copolymer  system  (31,  32).  When  the  solution  is  heated  (xabN  decreases),  the  phase  sequence 
can  be  detennined  by  a  diagonal  trajectory  towards  structures  dictated  by  the  block  copolymer  A 
fraction  since  the  solvent  becomes  less  selective  (30,  33).  This  diagonal  trajectory,  applied  to  the 
neat  copolymer  phase  diagram,  produces  the  correct  phase  sequence  upon  heating  for  the 
solution:  cylindrical  — >  spherical  — »  disorder  and  spherical  — >  disorder  at  intennediate  and  high 
solvent  dilutions,  respectively. 

The  morphology  of  the  ABA  +  B  blend  is  influenced  by  the  length  of  the  midblock  selective 
solvent.  Up  to  now  the  solvent  has  been  modeled  as  a  single  B  bead,  but  as  the  molecular  weight 
of  the  solvent  increases,  additional  beads  can  be  used  to  represent  the  solvent.  As  the  number  of 
beads  increase,  the  solvent  becomes  a  homopolymer  and  chain  stretching  becomes  an  important 
contribution  in  detennining  the  equilibrium  morphology.  Due  to  this  extra  interplay,  there  will  a 
transition  between  microphase  and  macrophase  separation  with  increasing  solvent  molecular, 
which  we  observed  for  our  gel  system.  We  considered  a  50/50  blend  (by  volume)  of  the  A3B15A3 
triblock  and  B  homopolymer-solvent,  since  the  spherical  morphology  was  verified  by  both  our 
model  and  experimentally  (29)  for  a  low  molecular  weight  solvent.  For  this  system,  we  studied 
the  phase  evolution  for  three  interaction  parameters,  XabN  =  46.9,  59.4,  73.5,  as  a  function  of 
increasing  homopolymer-solvent  bead  length  or  equivalently  molecular  weight.  At  a  low 
interaction  parameters  of  XabN  =  46.9,  where  our  A3B15A3  and  Bi  system  is  initially  disordered 
(figure  2a),  we  begin  to  observe  the  onset  of  phase  separation  to  a  spherical  phase  when  the 
length  of  the  homopolymer-solvent  is  increased  to  3  beads  (figure  2b)  and  to  an  ordered 
spherical  phase  at  5  beads  (figure  2c).  Macrophase  separation  is  observed  when  the 
homopolymer-solvent  is  equal  to  7  beads  (figure  2d),  where  a  disordered  phase  composed  of 
mostly  the  homopolymer-solvent  coexists  with  a  cylindrical  microphase.  Similar  behavior  is 
observed  for  larger  interaction  parameters  such  as  XabN  =  59.4  or  XabN  =  73.5,  where  the  A-B 
repulsion  is  substantial  enough  to  microphase  separate  the  original  A3B15A3  and  Bi  system  into 
the  spherical  morphology  (figure  3a).  At  these  two  interaction  parameters,  macrophase 
separation  occurs  when  the  homopolymer-solvent  is  represented  by  5  and  4  beads  at  XabN  =  59.4 
or  XabN  =  73.5,  respectively  (figure  3b).  Variations  in  the  number  of  beads  associated  with 
macrophase  separation  are  due  to  the  increased  contribution  of  interfacial  energy  to  the  overall 
free  energy  with  increasing  interaction  parameter.  Since  the  majority  of  the  homopolymer- 
solvent  is  expelled  from  the  triblock  copolymer  when  macrophase  separation  is  observed,  the 
triblock  microphase  separates  into  the  morphology  associated  with  the  neat  A3B15A3  triblock. 
This  is  not  observed  at  XabN  =  73.5,  where  the  neat  A3B15A3  triblock  takes  the  gyroid 
morphology.  This  discrepancy  may  be  due  to  packing  frustrations  associated  with  this  geometry, 
where  residual  homopolymer-solvent  in  the  microphase  can  stabilize  the  cylindrical  phase. 
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Figure  2.  Phase  evolution  with  increasing  B  homopolymer-solvent  length  for  a  50/50  blend  (by  volume)  of  the 

triblock  copolymer  and  homopolymer-solvent  for  !:.Ab/u  =  6.0  or  XabN  =  46.9.  Isosurfaces  are  shown  for 
the  A  component  (red  online),  where  the  B  component  is  excluded  for  clarity. 


Figure  3.  Phase  evolution  with  increasing  B  homopolymer-solvent  length  for  a 
50/50  blend  (by  volume)  of  the  triblock  copolymer  and  homopolymer- 
solvent  for  XabN  =  59.4  and  73.5,  where  macrophase  separation  is 
observed  when  the  solvent  is  represented  for  5  and  4  beads,  respectively. 

Isosurfaces  are  shown  for  the  A  component  (red  online),  where  the  B 
component  is  excluded  for  clarity. 

This  transition  has  been  extensively  studied  for  AB  diblock  copolymers,  where  the  same  general 
observations  can  be  applied  to  our  triblock  copolymer  system.  When  the  relative  length  of  the 
homopolymer-solvent  is  considerably  smaller  than  the  copolymer,  corresponding  to  the  solvent 
limit,  it  segregates  quite  evenly  between  the  two  phases  due  to  the  relatively  large  translational 
entropy  of  the  homopolymer-solvent,  where  the  partitioning  relieves  chain  stretching  of  both 
blocks  at  the  expense  of  increased  A-B  contacts.  Once  the  size  of  the  microdomain  is  smaller 
than  the  undisturbed  homopolymer-solvent,  constraining  it  within  the  microdomains  will  cause  a 
loss  of  conformational  entropy  due  to  the  compressed  configuration  of  the  homopolymer-solvent. 
Instead  the  homopolymer-solvent  demixes  from  the  copolymer  and  the  system  macrophase 
separates  into  a  coexisting  ordered  and  disordered  phase  formed  mostly  by  the  triblock 
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copolymer  and  homopolymer-solvent,  respectively  (28,  34-40).  Between  these  two  limits,  at 
intermediate  ratios,  the  translational  entropy  of  the  homopolymer-solvent  decreases  and  will 
segregate  into  the  microdomains  formed  by  the  same  species  in  order  to  reduce  the  interfacial 
energy  between  the  two  dissimilar  species. 

Since  the  B  end  blocks  are  constrained  to  an  interface,  the  mechanical  properties  of  triblock 
copolymers  are  different  from  that  of  the  analogous  AB  diblock.  Triblock  copolymers  can 
exhibit  a  looped  or  bridged  configuration,  where  the  B  interface  is  located  at  the  same  interface 
in  the  former,  while  in  the  latter  they  reside  on  different  interfaces  (16,  20,  41-43).  To  calculate 
elastic  properties  for  our  system,  we  used  a  finite-element-based  mechanical  homogenization 
method  for  linear  elasticity,  where  the  details  of  the  model  can  be  found  in  references  44  and  45. 
This  method  was  used  to  calculate  the  relative  values  of  the  Young’s  (E)  and  shear  (G)  moduli 
by  Andzelm  et  al.  (44)  for  a  SIBS  triblock  copolymer  melt  and  requires  two  input  parameters  for 
each  species:  the  Young’s  modulus  and  the  Poisson  ratio,  which  were  taken  to  be  3.34  (0.001) 
and  0.38  (0.499)  for  the  glassy  A  (rubbery  B)  phase.  The  h  nonn  of  the  Young’s  and  shear 
moduli  for  XabN  -73.5  as  a  function  of  ®  are  shown  in  figure  4  for  the  ordered  microphases  and 
were  calculated  from  the  components  of  the  moduli  ((Exx,  Eyy,  Ezz)  or  (Gxy,  Gxz,  Gyz))  in  the 
transverse  (x),  vertical  (y),  and  longitudinal  (z)  directions.  In  addition,  the  degree  of  monotropy, 
Ezz/Eyy  is  plotted  as  a  function  of  solvent  volume  fraction,  where  deviates  from  unity  indicate 
departure  from  isotropy.  Only  one  anisotropy  factor  is  needed  since  Exx  ~  Eyy  for  the 
morphologies  considered.  Figure  4.  The  h  nonn  of  the  (top)  Young’s  and  (middle)  shear  moduli 
as  well  as  our  (bottom)  anisotropy  factor  which  we  take  as  the  ratio  of  the  Young’s  modulus  in 
the  z-  and  x-directions  at  XabN  =  73.5  as  a  function  of  solvent  volume  fraction. 

indicates  that  the  spherical  (®s  =  0.5-0. 6)  morphology  is  isotropic  (Exx  -  Eyy  -  Ezz),  while 
monotropy  is  observed  in  the  cylindrical  (®s  =  0. 1-0.4)  and  gyroid  phase  (®s  =  0)  (Exx  ~  Eyy  < 
Ezz). Since  the  gyroid  phase  exhibits  cubic  symmetry,  there  should  be  no  directional  dependence 
in  the  Young’s  modulus  (46),  and  the  material  is  isotropic  rather  than  anisotropic  as  our 
calculation  suggest.  This  deviation  is  due  to  the  defective  nature  of  our  structure  and  our  shearing 
process  which  promotes  anisotropy.  In  tenns  of  the  mechanical  properties,  the  Young’s  modulus 
tends  to  decreases  with  increasing  solvent  concentration,  indicating  a  reduction  in  stiffness  and 
resistance  to  stretching  with  dilution.  The  shear  modulus  indicates  that  resistance  to  volume 
preserving  shear  decreases  in  the  following  phase  sequence:  gyroid  — >  spherical  — >  cylindrical. 
The  Young’s  and  shear  moduli  for  the  two  ordered  phases  shown  in  figure  2c  and  d,  were  also 
calculated  to  be  (0.283  GPa,  0.290  GPa,  0.288  GPa)  and  (0.093  GPa,  0.094  GPa,  0.097  GPa), 
respectively,  for  the  microphase  separated  spherical  phase  and  (0.241  GPa,  0.217  GPa,  0.492 
GPa)  and  (0.048  GPa,  0.047  GPa,  0.074  GPa),  respectively,  for  the  macrophase  separated 
cylindrical  phase.  Thus,  when  the  B  homopolymer-solvent  macrophase  separates  from  the 
triblock  copolymer,  the  system  becomes  anisotropic  in  the  z-direction,  where  the  system  exhibits 
less  resistance  to  shear  compared  to  the  spherical  phase. 
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Figure  4.  The  12  norm  of  the  (top)  Young’s  and  (middle)  shear 

moduli  as  well  as  our  (bottom)  anisotropy  factor  which  we 
take  as  the  ratio  of  the  Young’s  modulus  in  the  z-  and  .in¬ 
directions  at  XabN  =  73.5  as  a  function  of  solvent  volume 
fraction. 


4.  Conclusion 


In  conclusion,  we  investigated  the  phase  behavior  for  an  A3B15A3  triblock  immersed  in  a  Bi 
selective  solvent  as  a  function  of  solvent  volume  fraction  and  A-B  repulsion  in  order  to  model  an 
experimental  system  composed  of  poly[styrene-b-(ethylene-co-butylene)-b-styrene]  in  a 
midblock  selective  mineral  oil.  Solvent  swelling  of  the  microdomains  modifies  the  phase 
behavior  from  the  neat  triblock  copolymer,  where  the  gyroid  phase  becomes  unstable  with  small 
solvent  concentrations  and  the  spherical  morphology  becomes  the  sole  microphase  at  high 
dilutions.  The  effect  of  solvent  length  was  also  studied  for  a  50/50  (by  volume)  mixture  of 
copolymer  and  solvent,  where  a  transition  from  a  microphase  to  macrophase  separated  solvent 
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was  observed  after  a  critical  length  due  to  the  conformational  entropy  loss  associated  with 
confining  a  long  homopolymer-solvent  within  a  relatively  compressed  microdomain.  Finally, 
elastic  properties  were  calculated  using  a  finite-element-based  mechanical  homogenization 
method  for  the  gyroid,  cylindrical,  spherical,  and  macrophase  separated  cylindrical  phases. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


ARL 

U.S.  Army  Research  Laboratory 

C14 

hydrocarbon  oil  tetradecane 

DDFT 

dynamic  density  functional  theory 

DoD 

Department  of  Defense 

MS 

Material  Studios 

ODT 

order-disorder  transition 

SBS 

poly[styrene-butadiene-styrene] 

SCFT 

self-consistent  field  theory 

SEBS 

poly[styrene-(ethylene-co-butylene)-styrene] 

SIBS 

poly[styrene-isobutylene-styrene] 

SIS 

poly[styrene-isoprene-styrene] 

TPE 

thermoplastic  elastomer 
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